Hole burning in a nanomechanical resonator coupled to a Cooper pair box 
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Abstract 

We propose a scheme to create holes in the statistical distribution of excitations of a nanomechanical resonator. It 
employs a controllable coupling between this system and a Cooper pair box. The success probability and the fidelity 
are calculated and compared with those obtained in the atom-field system via distinct schemes. As an application we 
show how to use the hole-burning scheme to prepare (low excited) Fock states. 
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1. Introduction 

Nanomechanical resonators (NR) have been studied in a diversity of situations, as for weak force detections [1], 
precision measurements [2], quantum information processing [3], etc. The demonstration of the quantum nature 
of mechanical and micro mechanical devices is a pursued target; for example, manifestations of purely nonclassical 
behavior in a linear resonator should exhibit energy quantization, the appearance of Fock states, quantum limited 
position-momentum uncertainty, superposition and entangled states, etc. NR can now be fabricated with fundamental 
vibrational mode frequencies in the range MHz - GHz [4, 5, 6]. Advances in the development of micromechanical 
devices also raise the fundamental question of whether such systems that contain a macroscopic number of atoms 
will exhibit quantum behavior. Due to their sizes, quantum behavior in micromechanical systems will be strongly 
influenced by interactions with the environment and the existence of an experimentally accessible quantum regime 
will depend on the rate at which decoherence occurs [7, 8]. One crucial step in the study of nanomechanical systems 
is the engineering and detection of quantum effects of the mechanical modes. This can be achieved by connecting 
the resonators with solid-state electronic devices [9, 10, 11, 12, 13], such as a single-electron transistor. NR has 
also been used to study quantum nondemolition measurement [13, 14, 15, 16], quantum decoherence [12, 17], and 
macroscopic quantum coherence phenomena [18]. The fast advance in the tecnique of fabrication in nanotecnology 
implied great interest in the study of the NR system in view of its potential modern applications, as a sensor, largely 
used in various domains, as in biology, astronomy, quantum computation [19, 20], and more recently in quantum 
information [3, 21, 22, 23, 24, 25, 26] to implement the quantum qubit [22], multiqubit [27] and to explore cooling 
mechanisms [28, 29, 30, 31, 32, 33], transducer techniques [34, 35, 36], and generation of nonclassical states, as 
Fock [37], Schr6dinger-"cat" [12, 38, 39], squeezed states [40, 41, 42, 43, 44], including intermediate and other 
superposition states [45, 46]. In particular, NR coupled with superconducting charge qubits has been used to generate 
entangled states [12, 38, 47, 48]. In a previous paper Zhou and Mizel [43] proposed a scheme to create squeezed states 
in a NR coupled to Cooper pair box (CPB) qubit; in it the NR-CPB coupling is controllable. Such a control comes 
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from the change of external parameters and plays an important role in quantum computation, allowing us to set ON 
and OFF the interaction between systems on demand. 

Now, the storage of optical data and communications using basic processes belonging to the domain of the quan- 
tum physics have been a subject of growing interest in recent years [49]. Concerned with this interest, we present 
here a feasible experimental scheme to create holes in the statistical distribution of excitations of a coherent state pre- 
viously prepared in a NR. In this proposal the coupling between the NR and the CPB can be controlled continuously 
by tuning two external biasing fluxes. The motivation is inspired by early investigations on the production of new 
materials possessing holes in their fluorescent spectra [50] and also inspired by previous works of ours, in which we 
have used alternative systems and schemes to attain this goal [51, 52, 53]. The desired goal in producing holes with 
controlled positions in the number space is their possible application in quantum computation, quantum cryptography, 
and quantum communication. As argued in [52], these states are potential candidates for optical data storage, each 
hole being associated with some signal (say YES, |1), or |+)) and its absence being associated with an opposite signal 
(NO, |0), or |-)). Generation of such holes has been treated in the contexts of cavity-QED [53] and traveling waves 



2. Model hamiltonian for the CPB-NR system 

There exist in the literature a large number of devices using the SQUID-base, where the CPB charge qubit consists 
of two superconducting Josephson junctions in a loop. In the present model a CPB is coupled to a NR as shown 
in Fig. (1); the scheme is inspired in the works by Jie-Qiao Liao et al. [23] and Zhou et al. [43] where we have 
substituted each Josephson junction by two of them. This creates a new configuration including a third loop. A 
superconducting CPB charge qubit is adjusted via a voltage V\ at the system input and a capacitance C\ . We want the 
scheme ataining an efficient tunneling effect for the Josephson energy. In Fig.(l) we observe three loops: one great 
loop between two small ones. This makes it easier controlling the external parameters of the system since the control 
mechanism includes the input voltage V\ plus three external fluxes <&{€), <5(r) and <S> e (t). In this way one can induce 
small neighboring loops. The great loop contains the NR and its effective area in the center of the apparatus changes 
as the NR oscillates, which creates an external flux O e (f) that provides the CPB-NR coupling to the system. In this 
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Figure 1: Model for the CPB-NMR coupling. 



work we will assume the four Josephson junctions being identical, with the same Josephson energy 



Ej, the same being 
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assumed for the external fluxes <t>(f) and <t>(r), i.e., with same magnitude, but opposite sign: 0(f) = -<t>(r) = <$)(x). In 
this way, we can write the Hamiltonian describing the entire system as 



H = aid? ft. + 4E C |iVi - ijo% -4£°cos^^jcos|^jd-. t , (1) 

where a ' (a) is the creation (annihilation) operator for the excitation in the NR, corresponding with the frequency at 
and mass m; E a , and E c are respectively the energy of each Josephson junction and the charge energy of a single 
electron; C\ and stand for the input capacitance and the capacitance of each Josephson tunel, respectively <t>o = 
h/2e is the quantum flux and N\ — C\V\j2e is the charge number in the input with the input voltage V\. We have 
used the Pauli matrices to describe our system operators, where the states \g) and \e) (or and 1) represent the 
number of extra Cooper pairs in the superconduting island. We have: <x- = \g) (g\ - \e) (e\, & x = \g) (e\ - \e) (g\ and 
Z?c = e 2 /(C 1+ 4C5). 

The magnectic flux can be written as the sum of two terms, 

<D e = <Di + Btx , (2) 

where the first term <$>b is the induced flux, corresponding to the equilibrium position of the NR and the second term 
describes the contribution due to the vibration of the NR; B represents the magnectic field created in the loop. We 
have assumed the displacement x described as x = Xo(a^ + a), where xq = y/mZj/2 is the amplitude of the oscillation. 
Substituting the Eq.(2) in Eq.(l) and controlling the flux we can adjust cos f^H = to obtain 

H = utfa + 4E C (nj - i j & z - 4-E 1 } cos sin (^p j & x , (3) 

and making the approximation nB£x/<f>o « 1 we find 

H — wfi'fi + — o)q& z + Ao(a' + a)& x , (4) 

where the constant coupling A$ = -4£'5 cos {^~) ( ~ ) an d tn e effective energy a>o = SE C (Ni - . In the rotating 
wave approximation the above Hamiltonian results as 

H — (jja^a + — (jl)q& z + Aq(& + ci + fl'<r_). (5) 

Now, in the interaction picture the Hamiltonian is written as, Hj = U^HUQ-ihU^^-, where Uq = exp \—i (ata'a + ^y^) 
is the evoluion operator. Assuming the system operating under the resonant condition, i.e., a> = coq, and setting 
fr z = &+&- - and fr ± = (& x + i&yj /2 , with & y = (\e) (g\ - \e) (g\)/i the interaction Hamiltonian is led to the 

abbreviated form, 

Hi = J 8(a t d-_ + ad- + ), (6) 

where B = -Aq, &+ (d"_) is the raising (lowering) operator for the CPB. 

We note that the coupling constant B can be controlled through the flux (D x , which influences the mentioned small 
loops in the left and right places. Furthermore, we can control the gate charge A^i via the gate voltage V\ syntonized 
to the coupling. It should be mentioned that the energy a>o depends on the induced flux O v . So, when we syntonize 
the induced flux O t the energy changes. To avoid unnecessary transitions during these changes, we assume the 
changes in the flux being slow enough to obey the adiabatic condition. 

Next we show how to make holes in the statistical distribution of excitations in the NR. We start from the CPB 
initially prepared in its ground state \CPB) — \g) , and the NR initially prepared in the coherent state, \NR) = \a) .Then 
the state ^ that describes the intire system (CPB plus NR) evolves as follows 

\V NC (i)) = G(i)\g)\a), (7) 
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where U(t) = exp(-itHj) is the (unitary) evolution operator and Hi is the interaction Hamiltonian, given in Eq. (6). 
Setting &+ = \g) (e\ and <x_ = \e) {g\ we obtain after some algebra, 



U{t) = cos(/3fVa'a + Y)\g){g\ + cos(J3t ^f¥S) \e) (e\ 

, sinQ3fVfl t a+ 1) ., . sin^n/ga) 

-i — fl|g><e| -! = — fl'|e>(g|. (8) 

Vfl ' fi + 1 Vfl ' fl 

In this way, the evolved state in Eq.(7) becomes 

oo „ 

I^jvc(0> = c 2 >, -=[cos(w„T)|g,n) - /sin(w„r)|e,« + 1)], (9) 
£o Vn! 

where a>„ = /3 V« + 1, If we detect the CPB in the state \g) after a convenient time interval t\ then the state |TVc(0) 
reads 

oo n 

I^c(ti)> = m Y-=cos(co n T l )\n), (10) 
Vn! 

where 771 is a normalization factor. If we choose t\ in a way that /? V' J i + l T i = 7r /2, the component |nj) in the Eq.(10) 
is eliminated. 

In a second step, supose that this first CPB is rapidly substituted by another one, also in the initial state \g), that 
interacts with the NR after the above detection. For the second CPB the initial state of the NR is the state given in 
Eq.(10), produced by the detection of the first CPB in \g). As result, the new CPB-NR system evolves to the state 

00 n 

— =[cos(w„T 2 )cos(w„ri)|g,«> - /cos(w„Ti)sin(w„T 2 )|e,n + 1)]. (11) 

Next, the detection of the second CPB again in the state \g) leads the entire system collapsing to the state 

— — [cos(w„t 2 ) cos(w„ti) |n>], (12) 

where 772 is a normalization factor. In this way, the choice /3 V«2 + 1^2 = makes a second hole, now in the 
component |«2). 

By repeating this procedure M times we obtain the generalized result for the M — th CPB detection as 

00 Q n M 

\¥nc(tm)} = m ^ ~p ]~1 cos(co„Tj) \n) , (13) 
n=o V«! j =1 

where tj is the j - th CPB-NR interaction time. According to the Eq (13) the number of CPB being detected coincides 
with the number of holes produced in the statistical distribution. In fact the Eq (13) allows one to find the expression 
for the statistical distribution, P„ = |(«| v I , aic(tm)| 2 ; a little algebra furnishes 

(a 2 7«0nfi«»W;) 
p - I— . (14) 

" z: = o(« 2ffl /'«!)n^ic ^(w„T,)' 

To illustrate results we have plotted the Fig. (2) showing the controlled production of holes in the photon number 
distribution. 

The success probability to produce the desired state is given by 

00 M 

Ps = e-^ 2 ^(a 2m /m!) \\ cos 2 (aj n Tj). (15) 
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Figure 2: Holes in the photon number distribution, for a = 2.0, (a) at n\ = 4, for the 1" step; (b) at n\ =4 and «2 = 1, for the 2" d step; (c) at 
n\ = 4, m = 1 and n$ = 7, for the 3" 1 step. 



Note that the holes exhibited in Fig.(2)(a), 2(b), and 2(c) occur with success probability of 9%, 4%, and 0.3%, 
respectively. 

We can take advantage of the this procedure applying it to the engineering of nonclassical states, e.g., to prepare 
Fock states [60] and their superpositions [61]. To this end, we present two strategies: in the first we eliminates the 
components on the left and right sides of a desired Fock state \N), namely: \N - 1), \N-2), ...and \N + 1), \N + 2), in 
the second one, we only eliminate the left side components of a desired Fock state \N). In both cases, it is convenient 
to consider the final state of the NR as, 

oo n M 

|1Vc(tm)>' = i M Yj ^(-0 M ]~[ ^n+jTj) \n + M) , (16) 

which is easily obtained by detecting the Cooper pair box in the state \e). The success probability P' s to produce a 
Fock state \N) reads 

oo M 

P's = e-^ 1 Yji^'Vml) Y] sm 2 (oj„ +j Tj). (17) 

m=0 7=1 

In the first strategy, we prepare Fock states \N) with N - M, i.e., the phonon-number coincides with the number 
of CPB detections M. The fidelity of these states is given by the phonon number distribution at Pm associated with 
the state ^nc^m))' , 

ll!',sin'( \f]/*,) 



M 



(18) 



2Xo(tf 2 '7«D nf=, sin 2 ( Vn+lffr B )' 

We note that, in this case the fidelity coincides with the - th component of the statistical distribution Pn. The 
Fig. (3) shows the phonon-number distribution exhibiting the creation of Fock state |3), |4), and |5); all with fidelity of 
99%, for an initial coherent state with a = 0.6. 

In the second strategy, we prepare Fock states \N) with = 2M or 2M - 1 . The associated fidelity is also given 
by the Eq.(18). The Fig. (4) shows the phonon-number distribution exhibiting the creation of Fock states |3), |4), and 
|5), all them with same fidelity 99%, for an initial coherent state with a = 0.6. 



3. Conclusion 

Concerning with the feasibility of the scheme, it is worth mentioning some experimental values of parameters and 
characteristics of our system: the maximum value of the coupling constant p„,„ x as 45MHz, with B as 0, 1 T, t — 30/vm, 
x Q = 500/m and E { ] = 5GHz, with cj = 200nMHz. [6, 23, 42, 43, 55, 56, 57, 58, 59]. The expression choosing the 
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Figure 3: Phonon number distribution exhibiting the creation of Fock state: (a) |3> (P' s = 17%), (b) |4> (P' s = 11%), and (c) |5> (P' s = 7%); all with 
fidelity of 99% and initial coherent state with a = 0.6. 
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Figure 4: Phonon number distribution exhibiting the creation of Fock state: (a) |2> (P' s = 17%), (b) |3) (P' s = 1%), and (c) |4> (P' s = 0.3%); all with 
fidelity of 98% and initial coherent state with a = 0.6. 



time spent to make a hole, fi^rij + Itj = n/2, funishes tj ^ 0.3 ns, when assuming all the CPB previously prepared 
at t — 0. On the other hand, the decoherence times of the CPB and the NR are respectively 500 ns and 160 fis 
[58]. Accordingly, one may create about 1600 holes before the destructive action of decoherence. However, when 
considering the success probability to detect all CPB in the state \g), a more realistic estimation drastically reduces 
the number of holes. A similar situation occurs in [51, 52, 53], using atom-field system to make holes in the statistical 
distribution P n of a field state; in this case, about l//s is spent to create a hole whereas 1ms is the decoherence time 
of a field state inside the cavity. So, comparing both scenarios the present system is about 60% more efficient in 
comparison with that using the atom-field system. Concernig with the generation of a Fock state |A0, it is convenient 
starting with a low excited initial (coherent) state, which involves a low number of Fock components to be deleted via 
our hole burning procedure. According to the Eq. (18) when one must delete many components of the initial state to 
achieve the state |Af) this drastically reduces the success probability. As consequence, this method will work only for 
small values of N ( N < 5 ). 
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